Enhancing photoassociation rates by non-resonant light control of shape resonances 
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Photoassociation, assembling molecules from atoms using laser light, is limited by the low density 
of atom pairs at sufficiently short interatomic separations. Here, we show that non-resonant light 
with intensities of the order of lO^^W/cm^ modifies the thermal cloud of atoms, enhancing the 
Boltzmann weight of shape resonances and pushing scattering states below the dissociation limit. 
This leads to an enhancement of photoassociation rates by several orders of magnitude and opens the 
way to significantly larger numbers of ground state molecules in a thermal ensemble than achieved 
so far. 



I. INTRODUCTION 

The formation of ultracold molecules (T < 100 /iK) has 
experienced rapid progress over recent years [1]. At ul- 
tralow temperatures, extremely precise control over the 
molecules' dynamics can be exerted. This allows for 
transferring molecules into a single quantum state and 
even for controlling their reactivity [2]. Their rich inter- 
nal structure makes them sensitive probes in precision 
measurements of fundamental constants [3]. Possibly 
strong dipolar interactions facilitate their use in quan- 
tum computation and quantum simulation [4]. 

To produce molecular samples at ultralow temper- 
atures, molecules are typically assembled from ultra- 
cold atoms using magnetic-field controlled Feshbach reso- 
nances or laser light. As demonstrated in a series of spec- 
tacular experiments, ultracold molecules in their inter- 
nal ground state are obtained by subsequently applying 
Stimulated Raman Adiabatic Passage (STIRAP) [5-7] or 
broadband vibrational cooling [8]. However, in all these 
experiments, the number of ground state molecules that 
can be produced has so far been limited to about 10'*. 
When the molecules are formed by magnetoassociation, 
this small number is explained by the very low tempera- 
ture, of the order of hundreds of nano-Kelvin, and corre- 
sponding high phase space density, attained by e.g. evap- 
orative cooling. In the case of photoassociation, usually 
carried out in magneto-optical traps holding up to 10^*^ 
atoms at temperatures between 1 — 100 and a phase 
space density many orders of magnitude lower, only a 
small fraction of atom pairs resides at sufficiently short 
interatomic separations to be photoassociated. The num- 
ber of ground state molecules produced directly by pho- 
toassociation has thus also been limited to about 10'* [9]. 
A significant increase of the photoassociation rate is the 
problem that we address here. 

To this end, we suggest to control shape resonances 
using non-resonant light to enhance photoassociation 
rates. Non-resonant light couples to the anisotropic 
polarizability of an atom pair, shifting the position of 
shape resonances to lower energies [13] and increasing 
the resonance's thermal weight in an ultracold trap. Be- 



low we demonstrate an enhancement by three orders of 
magnitude in the photoassociation rate of strontium, a 
molecule which is currently attracting considerable atten- 
tion [11, 12, 14, 15]. Photoassociation relies only on the 
presence of optical transitions which usually are abun- 
dant. Shape resonances are ubiquitous for diatomics. 
The light-matter coupling via the anisotropic polarizabil- 
ity is of universal character, independent of a particular 
energy level structure, frequency of the light, or presence 
of a permanent dipole moment. Therefore our findings 
are easily extended to molecules other than Sr2. 

Our approach is related to Feshbach-optimized pho- 
toassociation [16] which also employs a scattering res- 
onance to increase the number of atom pairs that are 
quasi-trapped at short interatomic separations. Instead 
of tuning an existing resonance, one could also use reso- 
nances induced by an external field, for example, electric- 
field induced resonances for polar molecules [17, 18]. 
However, the electric fields required to achieve a signifi- 
cant enhancement (of about 1.3MV/cm in the most fa- 
vorable case) are currently experimentally unfeasible [19]. 
Unlike Feshbach-optimized photoassociation, shape res- 
onance control does not require the presence of a hy- 
perfine manifold of the atoms. Shape resonances occur 
naturally for partial waves with I > when a scatter- 
ing state becomes trapped behind the centrifugal barrier. 
Since shape resonances are quasi-bound states, they have 
very favorable free-to-bound overlaps required for effi- 
cient photoassociation; cf. the blue line in Fig. 1 (left). 
In some fortitious cases such as Rb2 [20, 21] or Cs2, the 
energy of a shape resonance almost matches the trap tem- 
perature, and large photoassociation rates are observed. 
Typically, however, the lowest energies at which shape 
resonances occur are far from the trap temperature, and 
their thermal weight is very small. In this case, pho- 
toassociation is limited by the atom pair density p(i?) 
at interatomic separations R < 200 ao [22]; i.e., it ad- 
dresses mostly s-wave atom pairs for which free-to-bound 
Franck-Condon factors are very small, cf. the grey line in 
Fig. 1 (left). This can be changed by applying a strong 
non-resonant field prior to and during photoassociation. 
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FIG. 1: (Color online) left: The photoassociation efficiency 
is determined by the free-bound Franck-Condon overlap be- 
tween the scattering states and weakly bound excited state 
levels (purple line), shown here for photoassociation of '**Sr2 
molecules near an intercombination line transition [10-12]. 
For s-waves (grey line, scaled by a factor of 10® compared 
to the blue line), the probability density at short internu- 
clear distances is very small, leading to low photoassociation 
efficiency. A shape resonance (blue line corresponding to a 
field-free rotational quantum number / = 8) is a quasi-bound 
state trapped behind the rotational barrier. Typically, its en- 
ergy, EsR, is too high for the shape resonance to carry any 
significant thermal weight in an ultracold cloud. Here we sug- 
gest to employ a strong non-resonant field (red arrow) which 
leads to rotational hybridization of the scattering states, ef- 
fectively moving the position of the shape resonance to lower 
collision energies, Ecou, and increasing the thermal weight of 
the shape resonance. 



II. THEORETICAL FRAMEWORK 

We assume the non-resonant light, with intensity /, 
and a photoassociation laser, with intensity I pa, to be 
linearly polarized along the z-axis of the laboratory fixed 
frame. Then, in the body fixed frame, the interaction of 
a pair of atoms with these two fields is described by the 
Hamiltonian 
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and j = g,e for the electronic ground and excited state, 
respectively. A.pA = h(ujpA — i-^at) denotes the detun- 
ing of the photoassociation laser from the atomic res- 
onance, eo and c are the electric vacuum permittivity 
and speed of light. The photoassociation laser inter- 
acts with the transition dipole moment, D{R), coupling 
the electronic ground state to an electronically excited 



state. Vg{R) and T4(i?) denote the respective poten- 
tial energy curves and Tp and L^/2/zi?^ are the vibra- 
tional and rotational kinetic energies. The non-resonant 
light intensity couples to the polarizability anisotropy, 
Aa^{R) = aji (i?) — a-'j_{R) {j = g, e), where ceyj_{R) are 
the parallel and perpendicular molecular polarizability 
components. The behaviour at large internuclear dis- 
tances is well approximated by Silberstein's formula in 
terms of atomic polarizabilities [23, 24]. We consider 
^®Sr atoms trapped at a temperature of 1 ^K. The po- 
tential energy curve for the electronic ground state was 
taken from Rcf. [25] and adjusted to yield a scatter- 
ing length of as = — 2ao [35]. The potential energy 
curves for the electronically excited A^E+, c'^n„, and 
a'^S+ manifold, spin-orbit coupling functions and transi- 
tion dipole matrix elements are found in Ref. [14]. We 
investigate photoassociation into the lowest level below 
the P asymptote that was previously observed 

{v' = —6) [10]. This level is well represented by an adia- 
batic potential energy curve 0,^ neglecting resonant spin- 
orbit coupling [15]. Its spontaneous emission lifetime was 
calculated to be 7„'=_6 = 7.6 ^s [15]. The atomic polar- 
izabilities for ground and excited state are taken from 
Refs. [26, 27]. The Hamiltonian (1) is represented by a 
mapped Fourier grid for the radial part and a basis set 
expansion in terms of Legendre polynomials for the an- 
gular part, taking advantage of the magnetic quantum 
number m being conserved. About A'^^ sa 2000 radial 
grid points and Imax ~ 20 are required to converge the 
calculations for T = f and / < 1.5 x 10^"W/cm^. 

Diagonalizing the Hamiltonian (1) for IpA = yields 
the bound levels and continuum states of the ground 
and excited electronic states in the presence of the non- 
resonant field. For 1^0, only the magnetic quantum 
number is conserved. The non-resonant field mixes dif- 
ferent partial waves of the same parity and I and /' are 
not good quantum number any more. For the sake of 
simplicity, we will label the field-dressed states by the 
field-free quantum numbers ?,m, and v',l',m', adding a 
tilde to I, v' and V to indicate that they are labels not 
quantum numbers. Note that for bound states, the field- 
dressed levels v' , V are adiabatically connected to the 
field-free quantum numbers w', /' even for very large in- 
tensities. The hybridization of the angular motion of the 
field-dressed wave functions is analyzed in terms of the 
rotational weights of the field-free partial waves. 
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where ^Pl^{R,9; Esp) = {R,0\i^?^^^{Esp)), and 

tl'f q{R,9; E) are the field-free scattering wavefunctions 
for partial wave I. The photoassociation rate coefficient 
in the presence of the non-resonant field is determined 
from the excited state bound levels and ground state 
continuum states, IV'^/r/ /) ^^'^ iV'f (E)) with {tpf (E)) 
normalized with respect to energy [29]. 

For an ensemble of atom pairs with a Maxwell- 
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Boltzmann velocity distribution, employing an isolated effcient is given by [28] 
photoassociation resonance, the photoassociation rate co- 



with "f.^jiiimi the spontaneous decay rate of level 
v'l'm', approximated by 7„'=_6,;'=o,m'=o, 7s = 
js{IpA, v' , V , m' , E, I, m, /) the stimulated emission rate, 
and 7 the total width, 7 = 7{;/p„, + 7s- Qt = 
(2TTfikBT/h^)^^^ denotes the translational partition func- 
tion, and A-,j,^, the detuning, A^,j,^, = E^,j,^, -hvpA- 
Since ®^Sr2 is a bosonic molecule with zero nuclear spin, 
the sum in Eq. (3) runs solely over even partial waves 
I in the electronic ground state. We consider only 
771 = m' = in Eq. (3), for which the effect of the 
non-rcsonant field is largest. The actual photoassocia- 
tion rate will be slightly higher than predicted here due 
to the neglected contributions from terms with m 7^ to 
the sum of Eq. (3). For weak photoassociation intensi- 
ties, the stimulated emission rate can be approximated 
by Fermi's golden rule, 

js{IpA,v\ l\ m', E, /, m, /) w 

^^IW'^,,„,l^(i?)cos^|^;^^(£;))P, 

i.e., by an integral over R and 0, ensuring the correct 
selection rules. 



III. RESULTS: PHOTOASSOCIATION OF SR2 

The photoassociation rate, Kj,i{IpAtOJpAtI^T)^ is 
shown in Fig. 2 for v' = —6 and weak (top) and strong 
(bottom) non-resonant fields. The black solid line shows 
the typical rotational progression of a photoassociation 
spectrum for / = 0. The largest peak is observed for 
Z' = 1, reflecting the pure s-wave character of the ther- 
mal cloud at 1 /iK. A weak non-resonant field (Fig. 2, 
top) shifts the photoassociation peaks. The largest shift 
occurs for V = 1. For the optical lattice of Ref. [12] with 
/ ~ 10''W/cm2, the shift amounts to about 3 MHz. An 
increase of the peak heights is only observed for T = 3 
and I' = 5. It is rationalized in terms of rotational hy- 
bridization in the excited state levels. Generally, bound 
levels are more strongly affected by the non-resonant field 
than scattering states. The I' = 3 and T = 5 states have 
thus some contribution from the I' ^ 1 partial wave which 
renders photoassociation from s-wave atom pairs into the 
field-dressed I' = 3 and /' = 5 bound state allowed, en- 
hancing the corresponding photoassociation peak. While 
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FIG. 2: (Color online) photoassociation lines (absorption rate 
coefficient) for photoassociation into v' = —6 for weak (top) 
and strong (bottom) non-resonant fields /. The intensity of 
the photoassociation laser is set to IpA ~ 1 W/cm'^, and the 
trap temperature to T = 1 /iK. 



this increase amounts to several orders of magnitude, the 
total photoassociation rate for I' ~ 3 and T = 5 is still 
smaller than the I' = 1 photoassociation rate without 
non-resonant field. 

At high non- resonant light intensities (Fig. 2, bottom), 
the rotational progression observed for / = is replaced 
by the pendular progression. More remarkably, a large in- 
crease of the I' = 1 peak, amounting to about three orders 
of magnitude compared to the field-free I' — 1 photoas- 
sociation rate, is observed in Fig. 2 for non-resonant field 
intensities of the order of 10^°W/cm^. For T ~ 1 /xK, 
the I = 12 shape resonance acquires its highest thermal 
weight close to / ~ 10^°W/cm^, cf. Fig. 3 showing the 
resonance position as a function of non-resonant field in- 
tensity. The remarkable enhancement of the I' = 1 pho- 
toassociation rate in Fig. 2 is thus related to the field- 
dressed I = 12 shape resonance. This is further ana- 
lyzed by inspection of the rotational weights, Eq. (2), of 
the field-dressed resonance wavefunction shown in Fig. 4. 
The lower (upper) panel of Fig. 4 displays the compo- 
nents for low (high) I. Note that the coupling mixes only 
partial waves of the same parity. For / = 0, the resonance 
wavefunctions are pure Z = 4, I = 8 and / = 12 states, 
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FIG. 3: Energy of the I = 4, 1 = 8, I = 12 and [ = 16 shape 
resonances vs non-resonant field intensity. 
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FIG. 4: (Color online) Contribution of different partial waves 
I to the resonance wavefunction for the ? = 4, 8, 12 shape res- 
onances as a function of the non-resonant intensity. The low-Z 
contributions are shown in the lower panel, the high-i contri- 
butions in the upper panel. For high non-resonant field in- 
tensities, all resonance wavefunctions acquire predominantly 
s-wave character. 

respectively. As the non-resonant field is switched on, a 
substantial amount of first I — 2 and then Z — 4 is mixed 
in to the field-free states. For large non-resonant field 
intensities, the trend of mixing in lower-Z components is 
observed all the way down to Z = (red diamonds). Thus, 
photoassociation from the field-dressed resonance wave 
function into V = 1 becomes possible for large /, ex- 
plaining the large enhancement of photoassociation into 
an excited state level with /' = 1 observed in the lower 
panel of Fig. 2. 

While the dependence on the non-resonant intensity 
in Fig. 4 is smooth for the / = 4 resonance, discontin- 
uous behavior is observed for ^ = 8, Z = 12. Each of 
the features in the apparantly complicated dependence 
for Z = 8, Z = 12 can be rationalized in terms of non- 
adiabatic behavior of interacting resonances: The first 
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FIG. 5: (Color online) I' = 1 peak in the presence of the field 
normalized with respect to the largest field-free peak (with 
I' = 1) in the rate coefficients for photoassociation into v' = 
—6 (T = 1 ^iK). Enhancement of the photoassociation rate of 
about three orders of magnitude is found when the position 
of a shape resonance comes close to the trap temperature. 



kink in the dependence of the I — 8 resonance is observed 
at 8 X 10^ W/cm^. It is due to the interaction with a Z ~ 6 
shape resonance [36] . A second discontinuous feature for 
Z = 8 is observed around 4 x lO^W/cm^. It is caused 
by the interaction of the / = 8 and I ~ 16 resonances 
close to their crossing, see Fig. 3. The most striking non- 
adiabatic features occur for both I = 8 and Z = 12 reso- 
nances between 8 x lO^W/cm^ and 1 x 10^°W/cm^. It 
is due to the non-adiabatic interaction between these two 
resonances which cross near 8.5 x lO^W/cm^ and with 
other scattering states, introducing a strong mixing of 
the resonance wavefunctions. The non-adiabatic behav- 
ior of the shape resonances observed in Fig. 4 underlines 
the requirement to treat the fully coupled rotational- 
vibrational dynamics. 

The largest increase of the V = 1 peak in Fig. 2 is 
found for I = 1.10 x 10^°W/cm^ (red dotted curve). 
A detailed analysis of the increase as a function of the 
non-resonant light intensity is presented in Fig. 5 which 
shows the height of the I' = 1 peak normalized with re- 
spect to the largest field-free {V = 1) peak. Most impor- 
tantly. Fig. 5 demonstrates an enhancement of the pho- 
toassociation rate between two to three orders of mag- 
nitude over a broad range of non-rcsonant field inten- 
sities. The narrow peaks at / 1.01 x 10^"W/cm^ 
and / = 1.225 x IQiOW/cm^ are due to the [ = 12 
and I = 8 shape resonances, cf. Fig. 3. These peaks 
refiect how the shape resonances pass through the en- 
ergy window that corresponds to the atomic cloud's ther- 
mal width. The broad features are attributed to regu- 
lar scattering states which are also affected by the non- 
resonant field. For sufficiently large J, a scattering state 
is pushed below the dissociation threshold (without sym- 
metry restriction, an even parity state is followed by an 
odd parity state) [18]. Before this happens, the scat- 
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FIG. 6: (Color online) Energy levels (a) and projections 
(b) of the hybridized v" = -3,1" = (red) and v" = 
—4, 1" = (blue) ground state levels for adiabatic (a) or sud- 
den (b) transfer into their field-free counterparts. For a slow 
switch-off, non-adiabatic transitions are not expected and all 
molecules are obtained with I" — (a), while for a sudden 
switch-off, about 20% of the molecules have I" = and about 
50% I" = 2. 



tcring wave function is compressed to short internuclear 
distances, yielding an increase in the photoassociation 
rate. This is observed in Fig. 5 for intensities around 
I = 1.085 X 10-^°W/cm^. The corresponding state be- 
comes bound at / w 1.22 x 10^°W/cm^. Also the peak 
near / = 5 x 10^ W/cm^ is caused by a state being pushed 
below the dissociation threshold. However, the enhance- 
ment is much larger for / ~ 1.085 x 10^° W/cm^ since this 
state becomes bound inbetween two shape resonances. 

The sharp features of the peak enhancement shown in 
Fig. 5 clearly reflect the I = rotational weight of the 
I — 8 and / = 12 shape resonances. The small peak 
at 1.01 W/cm^ is connected to the |ci=oP weight of the 
I = 12 shape resonance; and the sharp rise and drop of 
the normalized peak rate in Fig. 5 at 1.225 W/cm^ re- 
flects the behavior of the |q=oP weight for Z = 8 shape 
resonance before it becomes bound. The I = 16 shape 
resonance, also shown in Fig. 3 suffers almost no hy- 
bridization of its angular motion, with |c;=iep remaining 
as high as 0.9 for / = 6.1 x lO^W/cm^, just before the 
resonance becomes bound. Therefore the dipole coupling 
to v' = —6, / = 1 remains rather small, and no enhance- 
ment due to ^ = 16 shape resonance is observed in Fig. 5. 

In absolute numbers, the maximum peak in Fig. 5 cor- 
responds to a rate of 3 x 10~^^cm'^/s, comparable to the 
photoassociation rates of Cs2 molecules in weakly bound 
levels of the 0~(P3/2) excited state [30]. The important 
difference to Cs2 is due to the predominantly long- 
range behavior of the electronically excited state. This 
allows almost all of the photoassociatcd molecules to de- 
cay into bound ground state levels and efficient subse- 
quent Raman transfer to low lying vibrational states. 

The bound states of both ground and excited elec- 
tronic state are characterized by strong alignment in the 



presence of the non-resonant field [31]. For example, 
(cos'^e) > 0.94 for I > 10i"W/cm2. In the presence 
of the non-resonant field, the photoassociated molecules 
spontaneously decay into hybridized levels of the elec- 
tronic ground state. The branching ratios are calcu- 
lated from the transition dipole moments [15] of the field- 
dressed rovibrational states. For / = 1.1 x 10^°W/cm^, 
v' = —6, 1' = 1 shows a branching ratio of 31% to 
i" = -3,1" = and 43% to v" -4, [" com- 
pared to field-free branching ratios of 51% and 41% for 
v" = -3,1" = and I" = 2, respectively [37]. Once 
molecules in the electronic ground state are produced, 
the non-resonant field has to be switched off. If this 
is done adiabatically, the hybridized ground state lev- 
els v" -3,1" ^ and v" = -4,1" are directly 
connected with v" = -3, 1" = and v" = -4, /" = 0. A 
slow switch-off of the non-resonant field transfers these 
hybridized levels to their field-free counterparts unless 
avoided level crossings occur. For v" = —3, 1" = and 
v" — —A, I" — 0, crossings may be found only for very 
highly excited rotational states from lower vibrational 
bands, I" > 20, for which the non-adiabatic couplings 
with I" = are small, cf. Fig. 6(a). A slow switch- 
off of the non-resonant field thus produces ground state 
molecules with v" = —3 or v" — —4 and I" = 0. If 
the molecules after spontaneous decay are subject to loss 
processes in the trap, it may not be possible to switch 
the field off adiabatically. A sudden switch-off projects 
the hybridized ground state levels on a superposition of 
to' = field-free rotational levels of the same vibrational 
band with about 50% /' 2, 20% V = 0, and 20% 
I' — 4 components for both v" — —3 and v" = —4, 
cf. Fig. 6(b). Thus, a large part of the molecules is 
transferred to ground vibrational levels with I" = for 
both slowly or suddenly switching off the field. The 
weakly bound ground state levels show large two-photon 
transition matrix elements for efficient Raman transfer 
to low-lying rovibrational levels of the electronic ground 
state [15]. Almost all of the molecules produced by pho- 
toassociation can thus be transferred to the rovibrational 
ground state, in contrast to all previous photoassociation 
schemes realized experimentally to date. 



IV. CONCLUSIONS 

We have shown that the photoassociation rate of Sr2 
molecules is enhanced by three orders of magnitude by 
applying an additional non-resonant field. This enhance- 
ment is due to shape resonances and scattering states 
becoming bound, causing a larger number of atom pairs 
to be quasi-trapped at sufficiently short interatomic sep- 
arations to be photoassociated. Since the photoassoci- 
ation rate is limited by the low pair density at or near 
the Condon radius [22, 32], applying strong non-resonant 
light during photoassociation overcomes the main obsta- 
cle toward forming larger numbers of molecules. The en- 
hancement of the photoassociation rate is accompanied 
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by strong hybridization of the angular motion, which we 
have fully accounted for in a rigorous treatment of the 
coupled rovibrational motion. 

Prior to photoassociation, the non-resonant field needs 
to be switched on slowly such that the position of shape 
resonances and scattering states is shifted adiabatically 
to lower energies. Then the photoassociation laser ad- 
dresses a thermal cloud that is modified by the non- 
resonant light, producing molecules in an electronically 
excited state. As a specific feature of photoassocia- 
tion near an intercombination line, the majority of pho- 
toassociated molecules decays spontaneously into bound 
ground state levels rather than re-dissociate [10, 15]. In 
the presence of the non-resonant field, the spontaneous 
emission involves hybridized levels. To obtain field-free 
ground state molecules, the non-resonant field needs to 
be switched off. Ideally, this is done adiabatically, yield- 
ing rotationless molecules in two ground state vibrational 
levels. Alternatively, a sudded switch-off leads to a super- 
position of ground state molecules with the three lowest 
rotational quantum numbers, 0, 2 and 4, occupied. 

The required non-resonant field intensities of the order 
of 10^" W/cm^ can be achieved by employing intra-cavity 
beams. Assuming a spot size of the order of 10 /zm, intra- 
cavity powers of up to 10* W arc needed. Such powers 
can realistically be produced, for example, by injecting 
100 W from a single mode telecom wavelength fiber laser 



into a cavity with a power build-up factor of one hundred. 

Enhanced photoassociation rates arc a crucial prereq- 
uisite to increase the number of ground state molecules 
that can be produced optically. For molecules other than 
Sr2, our scheme works best for heavy atoms with large 
polarizabilities and large scattering lengths. Promising 
candidates include, besides other even-isotope alkaline- 
earth dimers or Yb2, heteronuclear molecules comprised 
of one alkali atom and one alkaline-earth-like atom such 
as RbSr or RbYb [33]. For these species, where Feshbach 
association is either not available or likely not to be fea- 
sible [34] , non-resonant field control of shape resonances 
paves the way to an efficient all-optical production of 
ground state molecules. 
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